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Spin-orbit coupling plays an increasingly important role in the modern condensed matter physics. 
For instance, it gives birth to topological insulators and topological superconductors. Quantum 
simulation of spin-orbit coupling using ultracold Fermi gases will offer opportunities to study these 
new phenomena in a more controllable setting. Here we report the first experimental study of a spin- 
orbit coupled Fermi gas. We observe spin dephasing in spin dynamics and momentum distribution 
asymmetry in the equilibrium state as hallmarks of spin-orbit coupling. We also observe evidences 
of Lifshitz transition where the topology of Fermi surfaces change. This serves as an important first 
step toward finding Majorana fermions in this system. 



In the past decade, quantum simulations with ultracold 
atoms have investigated many fascinating quantum phe- 
nomena in a highly controllable and tunable way. Studies 
of ultracold Fermi gases with resonant interaction have 
shed lights on understanding strongly interacting fermion 
systems in nature, including neutron stars and quark- 
gluon plasma; simulating Fermi Hubbard model with 
optical lattices helps to understand the physical mech- 
anism of unconventional high-Tc superconductors. How- 
ever, until very recently, an important interaction has not 
been explored in cold atom systems, that is, the spin- 
orbit (SO) coupling. Recently using two-photon Raman 
process, SO coupled Bose-Einstein condensate has been 
realized in the laboratory which gives rise to new 
quantum phases such as stripe superfiuid 0, HI . In real 
materials, SO coupling plays an important role in many 
physical systems over a wide range of energy scale, from 
determining the nuclear structure inside a nuclei, and 
the electronic structure inside an atom, to giving birth 
to topological insulators in solid state materials j3, HI- 
Since all these systems are fermionic, from the viewpoint 
of quantum simulation it is desirable to experimentally 
realize SO coupled degenerate Fermi gases. The physical 
effects of SO coupling in a degenerate Fermi gas are quite 
different from those in a Bose system. In this work we 
report the first experimental study of a SO coupled de- 
generate Fermi gas. Evidences of SO coupling have been 
obtained from both the Raman induced quantum spin 
dynamics and the spin-resolved momentum distribution. 
With SO coupling, we also observe evidences for Lifshitz 
transition where the Fermi surface changes its topology 
as the density of fermion increases. This progress will en- 
able us to study stronger pairing and higher Tc enhanced 
by SO coupling in resonant interacting Fermi gases 
and topological insulator and topological superfiuid in a 
more flexible setup [lO|, luj in the near future. 

In ultracold atom systems, SO coupling is generated 
through atom-light interaction induced artificial gauge 
fields, and therefore it has two advantages. First, both 
the strength and the configuration of SO coupling are 



tunable via controlling the atom-light coupling; Secondly, 
the coefficient of SO coupling is naturally on the same or- 
der of the inverse of laser wavelength, and in a gaseous 
system it gives rise to a coupling strength as strong as 
the Fermi energy. In this regime SO coupling dramati- 
cally changes the density-of-state at Fermi energy and the 
topology of Fermi surface, which gives rise to many in- 
triguing phenomena in many-body systems [6-12], while 
such a regime is not easy to access in conventional solid 
state materials. 

In our experiment, a degenerate Fermi gas of 2 x 10^ 
^^K atoms in the lowest hyperfine state \F = 9/2, = 
9/2) state is first prepared in an optical dipole trap. 
The optical dipole trap is composed of two horizontal 
crossed beams of 1064 nm at 90^ along the x ±y direc- 
tion overlapped at the focus, as shown in Fig. HJa). The 
temperature of the Fermi gas is about 0.3-0.4 Tp (Tp 
the Fermi temperature) when the trap frequency reaches 
27r X (116, 116, 164) Hz along (x, z) direction (see Meth- 
ods for details). We can change the trap frequency adi- 
abatically to achieve different fermion density. A pair 
of Helmholtz coils (green ones in Fig. 1(a)) provides a 
homogeneous bias magnetic field along y (quantization 
axis), which is precisely controlled by a carefully designed 
scheme described in Ref. [iS^ to reduce the magnetic field 
drift and the magnetic noise. 

The method we used to generate SO coupling is the 
same as reported by the NIST group for the ^'^Rb Bose 
condensate [1]. In the ^^K system, two spin- 1/2 states are 
chosen as two magnetic sublevels |t) = \F = 9/2, mp = 
9/2) and |;) = \F = 9/2, mp = 7/2) of the F = 9/2 hy- 
perfine level. They are coupled by a pair of Raman beams 
with coupling strength Q. Two Raman lasers with the 
wavelength A = 773 nm and the frequency difference u 
counter-propagate along x axis and are linearly polar- 
ized along y and z directions, respectively, corresponding 
to TT and a polarization relative to quantization axis y 
(as shown in Fig. 1(a)). The recoil momentum = 
/co sin(6>/2), and recoil energy Er = k'^ /2m = h x 8.34 
kHz are taken as natural momentum and energy units. 
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where ko = 2h7r/X and = 180^ is the angle between 
two Raman beams. A Zeeman shift a;z/27r = 10.4 MHz 
between these two magnetic sublevels is produced by the 
homogeneous bias magnetic field at 31 G. When the Ra- 
man coupling is at resonance (at uj/27t = 10.4 MHz and 
two-photon Raman detuning S = fi{ujz — oj) ~ 0), the 
detuning between \F = 9/2, mp = 7/2) and other mag- 
netic sublevels like |F = 9/2, mp = 5/2) is about x 170 
kHz, which is one order of magnitude larger than the 
Fermi energy. Thus we can safely disregard other levels 
and treat this system as a spin- 1/2 system. Same as in 
the boson experiment, this scheme generates an effective 
single particle Hamiltonian as [l| 
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Here, p denotes the quasi- momentum of atoms, which 
relates to the real momentum k as k = p =F kr-Bx with ^ 
for spin-up and down, respectively. This Hamiltonian can 
be interpreted as an equal weight combination of Rashba- 
type and Dresselhaus-type SO coupling [ij. Finally, be- 
fore time-of- flight measurement, the Raman beams, opti- 
cal dipole trap and the homogeneous bias magnetic field 
are turned off abruptly at the same time, and a magnetic 
field gradient along y direction provided by loffe coil is 
turned on. Two spin states are separated along y direc- 
tion during the time-of- flight due to the Stern- Gerlach 
effect, and imaging of atoms along z direction after 12 
ms expansion gives the momentum distribution for each 
spin component. 

The main part of this paper is to study manifesta- 
tion SO coupling in a Fermi gas system. We first study 
the Rabi oscillation between the two spin states induced 
by the Raman coupling. This experiment, on one hand, 
measures the coupling strength and on the other hand, 
clearly manifests different effect of SO coupling compared 
to boson system. All atoms are initially prepared in the 
It) state. The homogeneous bias magnetic field is ramped 
to a certain value so that S = —4:E^^ that is, the k = 
state of It) component is at resonance with k = 2k^ex 
state of It) component, as shown in Fig. [IJb). Then 
we apply a Raman pulse to the system, and measure the 
spin population for different duration time of the Raman 
pulse. Similar experiment in boson system yields an un- 
damped and completely periodic oscillation, which can 
be well described by a sinusoidal function with frequency 

[M]. This is because for bosons, macroscopic number 
of atoms occupy the resonant k = mode, and therefore 
there is a single Rabi frequency determined by the Ra- 
man coupling only. While for fermions, due to the Pauli 
exclusion principle, atoms occupy different momentum 
states. Precisely due to the effect of SO coupling, the 
coupHng between the two spin states and the resulting 
energy splitting are momentum dependent, and atoms 
in different momentum states oscillate with different fre- 
quencies (as shown explicitly in Eq. ([2]) later). Hence, 
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FIG. 1: Experimental setup and Raman- induced quantum 
spin dynamics: (a) Schematic of experimental setup and the 
Raman coupling of two hyperfine levels of ^°K. (b) The en- 
ergy dispersion with S = —AEr. The system is initially pre- 
pared with all atoms in |9/2,9/2) state, (c) The population 
in 1 9/2, 7/2) as a function of duration time of Raman pulse. 
/cf = 1.9A:r and T/Tf = 0.30 for red circles, /cf = l.SS/cr and 
T/Tf = 0.35 for blue squares, /cf = 1.1/cr and T/Tf = 0.29 
for green triangles. The solid lines are theory curves with 
Q = 1.52^r- (d) Time-of-flight image (left) and integrated 
time-of-flight image (integrated along y) at different duration 
time for | t) (blue) and | \) (red). The parameters are the 
same as blue ones in (c) 



dephasing naturally occurs and the oscillation will be in- 
evitably damped after several oscillation periods. This 
process is very similar to the spin dynamics of a spin po- 
larized current ejected into a semiconductor. For semi- 
conductor spintronics, one would like to have a polarized 
spin current, however, because different electron has dif- 
ferent velocity and thus their spins precess in different 
way, the current will be unpolarized. Such spin dynam- 
ics has been extensively studied during recent decades 
[iBj . In our case, we can observe momentum-resolved 
spin dynamics with Stern- Gerlach technique and there- 
fore we can clearly reveal this physical process. In the 
Fig. [IJc), we show the momentum distribution for both 
spin components at several different moments. One can 
see the multiple-peaks feature in momentum distribution 
in Fig. md), which clearly shows the out-of-phase oscil- 
lation for different momentum states. 



Unlike bosons, the Raman coupling Q cannot be read 
out directly from the period of oscillation. In fact, the 
frequency depends on both atoms density and tempera- 
ture, and it is generally smaller than Q due to the effect 
of averaging over different momentum states. To deter- 
mine the value of Q from the measurements, we fix Ra- 
man coupling and vary atoms density by changing the 
total number of fermions or the trapping frequency, and 
we obtain several different oscillation curve, as shown in 
Fig. [IJd). Then we fit them to the theory with a single 
fitting parameter Q. Theoretically, for a non-interacting 
system, the population of |^) component is given by 

n;(k + 2fcre:,,r,t) =n|(k,r,0) J"^ ^ o ^ 

sin^ V(^x/^r/m)2 + 1^2/4 1, (2) 

where t is the duration time of Raman pulse, n^(k, r, 0) 
is the equilibrium distribution of the initial state in local 
density approximation, and temperature of initial cloud 
is determined by fitting the time-of-flight image to mo- 
mentum distribution of free fermions in a harmonic trap. 
From Eq. ([2]) one can see that the momentum distribu- 
tion along X direction of |^) component is always sym- 
metric respect to 2kr at any time, and the experimental 
data is indeed the case, as shown in Fig. [U^c). The the- 
oretical expectation of the total population in ||) com- 
ponent is given by Ni{t) = J (i^k(i^rn|(k, r, t), and in 
Fig. md), one can see there is an excellent agreement 
between the experiment data and theory, from which we 
determine Q = 1.52(5)£^r- Since our current experiment 
is performed in the weakly interacting regime with 5- wave 
scattering length = 169ao, we have verified that the 
interaction effect is negligible (see Method for details). 

Next, we focus on the case with ^ = 0, and study 
the momentum distribution in the equilibrium state. We 
first transfer half of ^^K atoms from ||) to |t) using radio 
frequency sweep within 100 ms. Then the Raman cou- 
pling strength is ramped up adiabatically in 100 ms from 
zero to its final value and the system is held for another 
50 ms before time-of-flight measurement. We have also 
varied the holding time and find the momentum distri- 
bution does not change, and thus we conclude that the 
system has reached equilibrium in the presence of SO 
coupling. Since SO coupling breaks spatial reflectional 
symmetry {x ^ —x and kx —kx)-, the momentum 
distribution for each spin component will be asymmet- 
ric, i.e. ncr(k) ^ ncr(— k), with a =t,i; While on the 
other hand, when 5 = {) the system still preserves time- 
reversal symmetry, which requires n^(k) = n|(— k). The 
asymmetry can be clearly seen in the spin-resolved time- 
of-fiight images and integrated distributions displayed in 
Fig. [2fa) and (b), where the fermion density is relatively 
low. While it becomes less significant when the fermion 
density becomes higher, as shown in Fig. Efc), because 
the strength of SO coupling is relatively weaker com- 
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FIG. 2: Momentum distribution asymmetry as a hallmark of 
SO coupling: (a-c) time-of-flight measurement of momentum 
distribution for both | t) (blue) and | \) (red). Solid lines are 
theory curves, (a) /cf = 0.9/cr and T/Tf = 0.8 (b) /cf = l.G/cr 
and T/Tf = 0.63; (c) /cf = l.S/cr and T/Tf = 0.57. (d-f): plot 
of integrated momentum distribution ncr(k) — ncr(— k) for the 
case of (a-c). 



pared to the Fermi energy. In Fig. [2fa-c) the integrated 
momentum distribution is fitted by the theoretical calcu- 
lation to determine the temperature and the chemical po- 
tential at the center of the trap. (See Method for detail). 
We find that the Raman lasers indeed cause additional 
heating to the cloud. Nevertheless, the temperature we 
find is within the range of 0.5 — O.STp [ifij, which is still 
below degenerate temperature. In Fig. [2jd-f), we also 
show Ucrikx) — Ti(j{—kx) to reveal the distribution asym- 
metry more clearly. These are smoking gun evidences of 
a degenerate Fermi gas with SO coupling. 

With SO coupling, the single particle spectra of Eq. [T] 
are dramatically changed from two parabolic dispersions 
into two helicity branches as shown in the inset of Fig. 
ID^a). Here, two different branches are eigenstates of "he- 
licity" s and the "helicity" operator describes whether 
spin (Tp is parallel or anti-parallel to the "effective Zee- 



man field" hp = (—17, 0, k^Px/"^^^) at each momentum, 
i.e. 5 = tTp -hp/ltTp -hpl. 5 = 1 for the upper branch and 
s = — 1 for the lower branch. The topology of Fermi sur- 
face exhibits two transitions as the atoms density varies. 
At sufficient low density, it contains two disjointed Fermi 
surfaces with s = —1, and they gradually merge into 
a single Fermi surface as the density increases to nd- 
Finally a new small Fermi surface appears at the cen- 
ter of large Fermi surface when density further increases 
and fermions begin to occupy 5 = 1 helicity branch at 
nc2. A theoretical ground state phase diagram for the 
uniform system is shown in Fig. [H^a), and an illustra- 
tion of the Fermi surfaces at different density are shown 
in Fig. IH^b). Across the phase boundaries, the system 
experiences Lifshitz transitions as density increases [l7| . 
which is a unique property in a Fermi gas due to Pauli 
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principle. At sufficiently low temperature, the derivative 
of the thermodynamic quantities like the compressibility 
will exhibit singularity in the critical regime around the 
transition point. 

Rigorously speaking Lifshitz transition only exists at 
zero temperature and at finite temperature it becomes 
a crossover. However, we can still obtain several consis- 
tent features that supports the existence of such a tran- 
sition at zero temperature. We fix the Raman coupling 
and vary the atoms density at the center of the trap by 
controlling total fermion number or trap frequency, as 
indicated by the red arrow in Fig. [3ja). The quasi- 
momentum distribution in the helicity bases can be ob- 
tained from a transformation of momentum distribution 
in spin bases as follows (See Method for the definition of 
Up and Vp): 

^+(P) = — 72 J (3) 



ul - vl 



n_(p) = 



vl-ul 



(4) 



In Fig. [3fcl-c5), we plot the quasi- momentum distribu- 
tion in the helicity bases for different atoms density. At 
the lowest density, the 5 = 1 helicity branch is nearly un- 
occupied, which is consistent with that the Fermi surface 
is below 5 = 1 helicity branch. The quasi-momentum dis- 
tribution of the s = —1 helicity branch exhibits clearly 
a double-peak structure, which reveals the emergence of 
two disjointed Fermi surfaces at s = — 1 helicity branch. 
As density increases, the double-peak feature gradually 
disappears, indicating the Fermi surface of s = —1 helic- 
ity branch finally becomes a single elongated one, as the 
top one in Fig. [3l^b). Here we define a quality of visibil- 
ity V = (nA — n^)/{nx + ^b), where nx is the density of 
s = —1 branch at the peak and ne is the density at the 
dip between two peaks. Theoretically one expects v ap- 
proaches unity at low density regime and approaches zero 
at high density regime. In Fig[3l^d) we show that our data 
decreases as density increases and agrees very well with 
a theoretical curve with fixed temperature T/Tp = 0.65. 
Moreover, across the phase boundary between SFS and 
DFS-1, there will be a significant increase of population 
on 5 = 1 helicity branch. In Fig. [3] (e), the fraction of 
atom number population at s = 1 helicity branch is plot- 
ted as a function of Fermi momentum /cp, which grows 
up rapidly nearby the critical point predicted in zero- 
temperature phase diagram. The blue solid line is a the- 
oretical calculation for Nj^/N with T/Tp = 0.65, and 
the small deviation between the data and this line is due 
to the temperature variation between different measure- 
ments. Both two features are consistent with a Lifishitz 
transition. Recently, topological change of Fermi surface 
and Lifshitz transition have also been studied in Fermi 
gas in honeycomb optical lattices, where single particle 
spectrum exhibits Dirac point behavior [18]. A more 



accurate experimental determination of Lifshitz transi- 
tion point requires more well control of temperature in 
presence of Raman lasers and the information of local 
equation-of- state. In the near future, we also plan to 
bring the system close to a Feshbach resonance where 
the s-wave interaction becomes strongly attractive, and 
we will further cool the system below the superfiuid tran- 
sition temperature. There we expect to find Major ana 
fermion modes at the phase boundaries when Fermi sur- 
face topology changes 



In summary, we have for the first time studied effects of 
spin-orbit coupling in a ultracold atomic Fermi gas. We 
study Raman induced spin dynamics and reveal clearly 
the physical process of spin dephasing due to SO cou- 
pling. We measure the momentum distribution in the 
equilibrium state and find the momentum distribution 
asymmetry due to the broken of the refiection symme- 
try caused by SO coupling. We also find evidences for 
the change of Fermi surface topology and Lifshitz transi- 
tion from the quasi-momentum distribution in the helic- 
ity bases. All these features are unique manifestations of 
SO coupling in fermionic system, which are absent in pre- 
viously realized SO coupled boson condensate [1]. This 
research opens up an avenue toward rich physics of SO 
coupled Fermi gases. 

Methods. 

Experimental setup. The experimental setup (see Fig. 
1) is the same as discussed in (22l424j. which employs a 
mixture of ^^Rb and ^^K atoms. In short, a mixture of 
1 X 10^ ^^Rb atoms at the spin state \F = 2, = 2) and 
4 X 10^ "^^K atoms at |F = 9/2, ttif = 9/2) are precooled 
to 1.5 fiK by radio- frequency evaporative cooling in a 
quadrupole-Ioffe configuration (QUIC) trap, and then are 
transported to the center of the glass cell [25] in favor of 
optical access. Both species are loaded into the optical 
dipole trap and further evaporatively cooled in the op- 
tical dipole trap by evaporative cooling. At the end of 
above process, the 2 x 10^ ^^K atoms reaches quantum 
degeneracy at ~ O.STp. We use a resonant laser beam 
for 0.03 ms to remove the Rb atoms without heating ^^K 
atoms. 

Two laser beams for generating Raman coupling are 
extracted from a Ti-sapphire laser operating at the wave- 
length of 773 nm with the narrow linewidth single- 
frequency and focused at the position of the atomic 
clouds with 1/e^ radii of 200 /im. Two Raman beams 
are frequency-shifted with -100 MHz and -110.4 MHz by 
two single-pass acousto-optic modulators (AOM) respec- 
tively, and then are coupled into two polarization main- 
taining single-mode fibers in order to increase stability of 
the beam pointing and obtain better beam-profile qual- 
ity. The frequency difference uj of two signal generators 
for two AOMs is phase-locked by a source locking CW 
microwave frequency counters. To enhance the intensity 
stability of the two Raman beams, a small fraction of 
the light is sent into a photodiode and the regulator is 
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FIG. 3: Topological change of Fermi surface and Lifshitz 
transition: (a) Theoretical phase diagram at T = 0. k% — 
7i(37r^n)"^/^. "SFS" means single Fermi surface. "DFS" means 
double Fermi surface, (b) Illustration of different topology of 
Fermi surfaces as Fermi energy increases. The single parti- 
cle energy dispersion is drawn for small Q, in which red and 
green lines are for s = — 1 and s — \ helicity branches, respec- 
tively. Dashed blue line is the chemical potential. Red and 
green surfaces are Fermi surfaces for s = — 1 helicity branch 
and s — \ helicity branch, respectively, (c) Quasi-momentum 
distribution in the helicity bases. Red and green points are 
distributions for s = — 1 and s — \ helicity branches, re- 
spectively. /CF = O.g/cr, T/Tf = 0.80 for (cl); = 1.2A;r, 
T/Tf = 0.69 for (c2); /cf = 1.4/cr, T/Tf = 0.61 for (c3); 
A;f = 1.6/cr, T/Tf = 0.63 for (c4); /cf = 1.8/cr, T/Tf = 0.57 
for (c5). All these points are marked on phase diagram in (a), 
(d) Visibility v — (jik — tib) / {tia + ^b) decreases as k^/kr in- 
creases (A and B points are marked in (cl)). (e) Atom number 
population in s = 1 helicity branch N-^/N increases as k^/kr 
increases increases. In both (d) and (e), the blue solid line 
is a theoretical curve with T/Tf = 0.65, and the background 
color indicates three different phases in the phase diagram. 



used for comparing the intensity measured by the pho- 
todiode to a set voltage value from the computer. The 
non-zero error signal is compensated by adjusting the 
radio-frequency power in the AOM in front of the optical 
fiber. 

Calculating Momentum Distribution in a Trap. The 
energy eigenstates for the single particle Hamiltonian in 
Eq. (1) are the dressed states usually denoted as helicity 
bases with s = ±1, with |+,p) = p) + '^pl p) 



and |- 
- 



5 p) = '^pI P) ~ '^pI i-,v)-> where the coefficients 

are only functions 
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of px with Vr = kr/m. The energy dispersion for each 
branch of these dressed states are given by ^±,p = (p^ + 
k^)/{2m) ± Vfe^r - V2)^ + ^74. 

The trapping potential V^(r) is taken into account by 
local density approximation. Even in presence of Raman 
coupling, what can be measured in the time-of-flight im- 
age with Stern- Gerlach technique are still momentum dis- 
tributions of original hyperfine spin states. Theoretically 
they are given by 



dh 



(27rfi)3 

^3. 



ii^/(^p,+;r) + ^^/(^p,_;r) 



n^(p - kj-ex) -- 

[^p/(^P,+; r) + ^p/(^p,- ; r) 

where f{E]r) is the Fermi-Dirac distribution function 
with a local chemical potential /i(r) = /i — V{t), and fi is 
determined by the equation for total particle number 



N= d^p nt(p)+n;(p) 



(5) 



By fitting the measured data to the theoretical 
curve of integrated momentum profile ni^^aiPx) = 
J dpydpz ncr{p) with <j =^ or ^, we obtain the temper- 
ature of atoms in the experiment. 

Interaction Induced Change of Rabi Frequency. We 
shall also compare this experiment to the clock experi- 
ment in Fermi gases [26] , where the coupling between two 
components is momentum independent if the light inten- 
sity is assumed to be uniform. In that case, all atoms 
undergo Rabi oscillation with exact same frequency, and 
therefore remain as identical particles during the process. 
Hence, the s-wave interaction plays no role except when 
the spatial inhomogeneity is taken into account [26]. In 
contrast, in our case different atoms oscillate with dif- 
ferent frequencies, thus they immediately become distin- 
guishable particles as the oscillation starts, and can in- 
teract with each other via s-wave collisions. Since our ex- 
periment is performed in the weakly interacting regime, 
we include the interaction effect with mean-field theory. 

For a non-interacting system, the equation of motion 
for a spin Sp with quasi-momentum p is given by 



Sp X hp 



(6) 



where hp = {—ft^O^krPx/'m + S) is the momentum de- 
pendent "effective magnetic field" . The frequency of the 
spin rotation is given by |hp| = \/ {pxkr/ra + + Q^, 
and for 5 = —4:Ey^ it is just the frequency of popu- 
lation oscillation as shown in Eq. (j2j). For a weakly- 
interacting system, the interaction can be approximated 
by i^int = ^S-S at mean- field level, where g = Airh'^as/m 
is the coupling constant, as is the 5- wave scattering 
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length, and S = J2p total spin. Thus, an addi- 

tional mean-field term appears in the equation of motion. 



(hp 



2^8) 



(7) 



where the solution of Sp must be determined self- 
consistently. We numerically solve this equation of mo- 
tion to determine the spin dynamics. Using the exper- 
imental parameters, we find the frequency shift is only 
a few percent of which is beyond the measurement 
resolution of current experiment. 
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